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DYNAMICS OF LIQUID PROPELLANT GUN LOADING SYSTEMS 



Annual Progress Report, FY75 



I. INTRODUCTION AND BACKGROUND 

A. LIQUID PROPELLANTS FOR ADVANCED NAVY GUN SYSTEMS 

With continued trends toward more sophisticated weapons platforms 
to meet the projected threats of a hostile environment, improved gun 
systems will be required. These systems of the future will have to 
have a total weight less than existing gun systems of the same cali- 
ber, a less complicated logistics train, improved capabilities, and 
they will have to be cost effective. Various analytical and experi- 
mental investigations conducted over a period of approximately thirty 
years have demonstrated that the potential advantages of liquid propel- 
lant guns (LPG's) are sufficient to warrant their serious considera- 
tion as candidate elements in advanced naval gun systems. 

The following remarks are intended to illustrate some of these 
potential advantages. By eliminating the cartridge case, both air- 
borne and surface platforms will be able to carry more shots per 
pound and/or a more varied assortment of weapons. At the present 
time, shipboard guns using semi-fixed ammunition require a separate 
magazine for powder in relatively close proximity to the gun. By 
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using a liquid propellant, the propellant tanks could be located 
anywhere in the ship without great difficulty, thus allowing the 
ship designer more flexibility. With a standardized propellant, 
different caliber guns, which now use both semi-fixed and fixed 
ammunition, could be supplied from the same tank, eliminating the 
separate spaces needed to conform with ordnance storage compatibility 
requirements. Redistribution of propellants onboard could be accom- 
plished by means of pumping rather than manual labor or an elaborate 
mechanical transfer system. Propellant retrograde materials, e.g., 
pallets, powder cases, and powder cans, would be eliminated. Under- 
way replenishment and rearming evolutions would consist of taking on 
projectiles and liquid, thus significantly reducing time, cost, and 
effort. Ship stability might be improved by the elimination of propel- 
lant weight in ready-service spaces. A ship sustaining battle damage 
could easily improve her stability by redistributing propellant or 
pumping it over the side. Increased firing rates may be attainable. 

It is because of this prospectus and improved technology in the field 
of monopropellants and bipropellants for gun applications that a re- 
search and development effort has again been initiated in the liquid 
propellant gun (LPG) area of ordnance development. 

B. REVIEW OF PREVIOUS DEVELOPMENTS 

In order to establish a point of departure for current develop- 
mental efforts in LPG technology, it is useful to briefly review the 
experiences documented by previous investigators. A review of LPG 
history was conducted, but the following remarks refer only to un- 
classified literature. 
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Interest in LPG systems in the United States was evident as early 
as 1945 [1] when the Army Office of the Chief of Ordnance conducted 
some preliminary calculations of the performance expected from liquid 
propellants. It was concluded that artillery performance could be 
greatly improved by using liquid propellants instead of solid pro- 
pellants. Initial experimental investigation was in the area of bi- 
propellants and feasibility demonstrations were conducted with a 
smooth bore, caliber .50 gun. These early tests resulted in surpris- 
ingly high muzzle velocities which seemed to spur increasing interest. 
It was discovered that to inject the propellant properly, an injection 
pressure in excess of chamber pressure was required [1]. To circum- 
vent this problem, pre-loaded cased rounds were used. Large variations 
in the ballistic behavior without repeatability were noted [2]. 

In 1948, investigation into the feasibility of regenerative injec- 
tion of bipropellants was begun [3]. The regenerative injection tech- 
nique utilizes a differential area piston to extract energy from the 
burning propellant to power the injection process. It was determined 
that regenerative injection was feasible. However, continued study 
through 1950 with caliber .50 guns emphasized the need for high pres- 
sure static and dynamic seal improvement. It was also noted that 
lower mixing intensities gave lower muzzle velocities and excessive 
round-to-round variations were encountered with dilute bipropellants 
[4]. On the positive side, liquid propellants produced less bore wear 
than the existing solid propellants. 

Concurrently , work was started on a 37-mm gun, based on the caliber 
.50 work [5, 6]. The 37-mm gun had regenerative injection with 
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nitrogen pressure to start the process. Feasibility of using the 
regenerative process in 37mm guns was proven; however, further dev- 
elopment in the area of seals was again indicated. 

The Navy became interested and began investigation of LPG systems 
in 1949 with an emphasis on improved performance for aircraft guns 
[7, 8]. Larger guns with higher velocities at the same size and 
weight were also desired. Naval Ammunition Depot, Seal Beach, Calif- 
ornia, started research and development of a 40mm monopropellant gun. 
The monopropellant developed seemed to satisfy all of the desired 
propellant criteria, however, for their initial gun tests, a bipropel- 
lant was used in a pre-loaded casing. At the same time, separate 
research was being done in the area of propellants themselves [9]. 

It should be noted that, as of 1950, the majority of the inter- 
est in liquid propellant guns was in the Army. They were mainly in- 
terested in LPG technology vice development, with a major interest in 
the end application to antiaircraft weapons. The high muzzle veloci- 
ties and high rates of fire anticipated could readily satisfy the 
requirements for this role. Another use contemplated for LPG 1 s was 
aircraft machine guns, because of the high rates of fire and the 
elimination of cartridge cases. 

With the commencement of hostilities in Korea, the Army's LPG 
interest shifted to the area of infantry and tank weapons. The hyper- 
velocities and greatly increased number of rounds per reload were 
extremely desirable characteristics for infantry weapons. A program 
specifically designed to employ liquid propellants in an existing 
90mm gun as a tank weapon was started in 1950 [10]. This retrofit 
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program was aimed at improving the armor penetration of Army guns in 
Korea. The concept was to fit into a standard cartridge case the 
components of a hypergolic bi propellant system necessary to achieve 
muzzle velocities in excess of a standard round. Hypergolic bipropel- 
lants are designed to ignite spontaneously upon mixing. For a liquid 
propellant gun, this implies that no initiator is needed for ignition. 
These hypergolic bipropellants require an injection or mixing process. 

In the event of accidental contact of fuel and oxidizer, a fire would 
inevitably occur. Separate fuel and oxidizer containers were placed 
in the cartridge case with a mechanism to effect mixing. After two 
years of work, the design objectives had not been reached; only a half 
charge of propellant could be used without over-pressurizing the gun. 

Actual small-caliber design studies of caseless guns were started 
in the early 1950's. The major emphasis in this field was on aircraft 
machine guns. A caliber .60 gun was designed in 1953 incorporating a 
regenerative injector monopiston to introduce both fuel and oxidizer 
in equal amounts into the chamber [11 J. The design rate of fire was 
2,000 rounds per minutes. However, the projectile and propellant feed 
systems interacted in such a way that projectile feed could inhibit 
propellant loading. It was recommended that a revolver type feed 
system be investigated. The reproducibil ity of ballistic character- 
istics, misfire procedures, overpressurization of the gun, propellant 
ignition, and ullage problems continually plagued these early efforts 
in liquid propellant gun research. 

Until 1955, LPG systems consisted primarily of two general classes, 
preloaded and injection guns. Bipropellants and monopropellants have 
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been employed in each system with varying degrees of success. Pre- 
loaded guns, though simple in design, were extremely sensitive to 
loading density and ignition variations. Injection guns, though 
relatively complex, offered a positive means of combustion control. 

In an attempt to utilize the advantages of the various systems pre- 
viously investigated, a project was started in 1955 to combine the 
two systems into a single monopropellant system approaching the sim- 
plicity of a preloaded gun, yet having the advantage of combustion 
control through injection [12]. A 40mm model was designed, featuring 
regenerative injection. A low-pressure static, high-pressure dynamic 
seal configuration was adopted to circumvent earlier problems encountered. 
Marginal performance was obtained with the system. 

At the same time, the Navy funded research in the area of monopropel- 
lant and bipropellant injector systems for a 30mm gun for aircraft appli- 
cations. This design incorporated a breech bleed feature to reduce ullage. 

A 40mm development program was started as a prototype for a 3 to 5 
inch rapid fire gun [13]. The 40mm was selected for testing purposes as 
a caseless gun. This effort was apparently terminated after the gun was 
damaged. 

A number of projects were in progress in the mid-50's which were 
aimed at producing a 90mm tank gun [14, 15, 16, 17]. These programs 
continued to strive for hypervelocities by employing regeneratively 
injected bipropellants. An impressive increase in gun-tube life was 
anticipated-, however, the injection system was somewhat large for tank 
application. At this point (1957), development aimed at major-caliber 
end items for Army use appears to have stopped, although some small -arms 
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studies continued. Among the reasons for this curtailment of research 
would appear to be the movement to a peace-time economy with the cessa- 
tion of hostilities in Korea and the growth of missile technology. 

In 1968, programs were started by the Navy which addressed the 
development of a safer aircraft machine gun with improved performance 
[18]. The use of non-hypergol ic , preloaded, caseless liquid bipropel- 
lants was the course of action selected. After encouraging initial 
results, the program was continued with emphasis placed on developing 
a 20mm burst-fire gun. These programs utilized red fuming nitric acid 
(RFNA) as a propellant ingredient. One of these programs investigated 
a rotating four-barreled, cam operated system designed for a rate of 
fire of 1000 rounds per minute. It employed a pressurized injector 
feed system having an injection period of 11.3 milliseconds [19]. A 
computer simulation indicated the system would perform satisfactorily ; 
however, a working model was not constructed. A similar program con- 
ducted simultaneously also employed a multi-barreled design [20]. A 
single barrel test gun was manufactured and fired in single shot tests. 
In the multi-barreled design, the projectile would be seated by propel- 
lant pressure. A nitrogen pressurized accumulator system for the prop- 
ellant was employed to maintain propellant pressure. For the inert 
tests, water was used as the propellant simulator. 

C. SCOPE AND OBJECTIVES OF THE PRESENT STUDY 

The overall goal of this investigation is to both analytically and 
experimentally simulate and study the behavior of a liquid under con- 
ditions analogous to those existing in a rapid-fire LPG feed system. 
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The results of the study are intended to be useful in the identifica- 
tion of upper and lower bounds for such LPG design and performance 
parameters as time-to-load , injection supply pressure, injection sys- 
tem configuration (line lengths, diameters, accelerating masses, and 
viscous and minor losses), ullage, charge-to-mass ratio, caliber size, 
and projectile mass. In addition, the investigation has been designed 
to detect and characterize "anomalous" behaviors such as chamber pres- 
sure oscillations during loading, existence and consequences of cavita- 
tion in the liquid slug, and various sources of delay in the loading 
system. It has been felt that laboratory observation of effects such 
as these will prove to be essential in understanding the behavior of 
LPG prototype systems and future operational devices. 

D. APPROACH AND PROGRESS 

The investigation has taken the classical approach of analysis com- 
plemented by experimental tests. The analytical model has been con- 
structed with a view towards simplicity while maintaining flexibility 
to refine the analysis as required to meet design requirements and to 
incorporate new information as LPG technology advances. The analytical 
model in its current form assumes one-dimensional unsteady flow of an 
incompressible propellant. This latter assumption is viewed as the 
most critical deficiency in the analysis and is a recommended area 
for follow-on efforts. It will be seen, however, that the analytical 
model in its present form is adequate for use in obtaining prelimin- 
ary design guidance. 

The experimental portion of the study has been primarily designed 
to support and validate (hopefully) the analytical model. Accordingly, 



11 



no special effort has been made to emulate existing or envisioned LPG 
designs. Off-the-shelf hardware has been used in the construction of 
the laboratory LPG simulator so that the performance of the device is 
somewhat less than might be obtained from an optimally designed system. 
The analytical model accurately predicts this "off-design" performance 
and points the way toward the areas of redesign that are likely to be 
most beneficial. This redesign is another area recommended for the 
expenditure of future efforts. 

Progress to date may be summarized as follows: 

1. An analytical model of LPG loading systems has been constructed 
and exercised. The model has been shown to adequately predict the 
overall performance of LPG loading systems (time-to-load) and the 
dominant factors to be considered in the design of LPG systems. 

2. An LPG laboratory simulator has been designed and built and 
has been shown to be a useful tool for the evaluation of analytical 
predictions and for the identification of unpredicted loading system 
behavior. 

3. The analytical model, verified by the experiments, has been 
used to provide design guidelines for the estimation of LPG loading 
system performance as a function of various design parameters. 

These developments are discussed in detail in subsequent sections 
of this report. The report is concluded with a summary of results and 
recommendations for future work in the subject area of investigation. 
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II. ANALYTICAL MODEL 



Analytical efforts have been aimed at developing a model for the 
behavior of a liquid as it is pumped from an injector, through a system 
of connecting lines and valves, into the breech. Although continuing 
analytical work will endeavor to treat this behavior in some detail (in- 
cluding, for example, propagation of pressure waves and local occurrences 
of cavitation), the complexity of the process has dictated a quasi-lumped 
parameter approach for initial studies. As shown in the sketch below, 
the fluid is considered to be made up of three slugs occupying the in- 
jector, connecting lines, and breech. The present dynamic analysis takes 
into account the transfer of fluid from one region to another in response 
to a step input in pressure at the injector head end. 
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In order to conform with current design concepts, including the design 
of the laboratory simulator, the system is modeled in the "projectile- 
ram mode" in which the breech is the same diameter (nominally) as the 
bore. The model may be easily extended to the "flow-thru" mode of 
chamber filling by (in the first approximation) simply neglecting the 
mass of the slug. 

A. CONSTITUENT RELATIONSHIPS 

1. Injector inertia : The loss in pressure across the moving 

injector mass (see sketch) may be obtained from the force balance: 



P r\ ' P 1 A i = m pV Vi + k sp x i + F p 




( 1 ) 



2. Losses in the liquid : In this region there must be a 

pressure drop to overcome both the inertia and drag of the liquid. 
Thus, in the injector 

(L. - x.) x . 2 

AP i = Ptq pf, -V- !- 4r- (2) 



D. 

i 



in the connecting 1 ine 

Lx 2 

AP C = pL l*l + pf £ ~T 
and in the chamber 



X X 

aP c = px s x s + pf c W~ "T 

c 



( 3 ) 
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In addition, minor losses can be taken into account by writing 

2 



AP 



pX < 

ml ~ ? A ^i ~2 
AP, 



1 



E k • 

Y mi 



( 4 ) 



where k . = 
mi 



px. 



For example, for the drop due to the sudden contraction from the 
injector to the connecting line, we have (i =1) 

/ • v 2 

x. 



AP : iij, ^ . IfLjj 2 r! 
Ah l 2 2 x l 2 2 x s 1 



and k, 



ml 



3 

2 



Likewise, for orifices occurring in the connecting line, we would 
obtain (i = 2), for each such occasion. 




Invoking the incompressible assumption the velocities can all 
be written in terms of x s . Thus, since the volume flow from the 
injector must equal the inflow to the chamber, 





A, x, = nA^ = A c i s 


(5) 


whence 


A i *i - "Vt * Vs 


(6) 


and 


A i x 1 = Vs 


(7) 
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Combining Eqs. (2) through (7), the total pressure drop across the 
liquid may be written. 




3 . Pressure drop from reservoir to slug: 
By simple manipulation, we have 




so that after applying Eqs. (5) through (7) to Eq. (1), combining 
with Eq. (8), and collecting terms, the total drop from reservoir 
to slug is 



P -p' 

-ZyS- = 1 - P s = k 1 + k 2 x (H 1 + HgX) + k 3 X 2 (H 3 + H 4 x) 



+ k 4 x + k 5 x + k g 




(9) 
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(the charge- to-mass ratio) 




4. Drop across the slug : 

Again, writing a force balance, 

P'A = m x + k. x + k x + F +P k A 
sc ss fss sss s be 



or P $ = k 2 x + k^ x + kg x + kg (10) 

where: 

k* x k x F P. 

I, - fs sm . ss sm „ _s , _b 

7 P r A c ’ K 8 " P r A c ’ R 9 ' P r A c + P r ‘ 

Equations (9) and (10) are the governing differential equations 
for the slug position and slug-face pressure as a function of time. 
In subsequent sections, we have developed solutions for these 
equations under the initial conditions 

x (0) = x(0) = 0 
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B. SOLUTIONS TO THE GOVERNING EQUATIONS 

Equations (9) and (10) have been combined and programmed for 
solution on the analog computer. This mode of solution has been 
chosen because of the ease of parameter variation made possible by 
the analog method. The analog computer solution is outlined in 
Section IV of this report. In this subsection, several closed form 
solutions are presented as they are especially useful for the spec- 
ial cases treated. 

1. Solution of the semi-linear case . Combining Eqs. (9) and 
(10), we obtain, upon elimination of P, 

k 2 x(l+H 1 +H 2 x) + k 3 x 2 (H 3 +H 4 x) +(k 4 +k 7 )x + (k g +k g )x = 1 - (k 1 +kg+k g ) (11) 

This expression is non-linear due to the fluid friction term 
(preceded by k 3 ) and the terms accounting for the shift of fluid 
from the injector to the chamber (terms preceded by H g and H 4 ). An 
important simplification can be made if the system design parameters 
are such that 

H4 H g 

1 * TTHp* 1 

and if k g + k g << 1 . 

The first set of inequalities requires that the difference between 
the velocity and acceleration of the fluid in the chamber and of that 
in the injector is small or that piston mass and connecting line ef- 
fects dominate the system. The second inequality neglects spring 
effects in the piston and slug behavior. Both of these conditions 
represent real (but not necessarily desirable) situations and, in 
fact, are close approximations to the conditions of the experiments 
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described in Section III of this report. With these simplifications 
we have, 

k g ( 1 + H-j )x + k 3 x 2 + (k^+k^)x = 1 - (k-j+kg+kg) 
x(0) = x(0) = 0 

This is a Riccati differential equation with the solution 



where 



x = A In c0$h B ( tH ° ) - Ct 
cosh Bt 

o 

x = AB tanh B(t+t Q )- C 
x = AB 2 sech 2 B(t+t Q ) 



A = 



k 2 (l+H 1 ) 
k 3 H 3 



/~q~ r _ k 4 +k 7 
2k 2 (l+H 1 ) ’ " 2k 3 H 3 



(13) 



q = (k 4 +k 7 ) 2 + 4 [1 - (k^kg+kg)] (k 3 H 3 ) 




tanh 



( k 4 +k 7 ) 

~7T 



The variation of pressure P with time may be obtained by substituting 
Eq. (13) into Eq. (9) or Eq. (10), but this rather laborious procedure 
has been omitted in deference to the slavishness of the computer. 

A most useful application of Eq. (13) is to use it to estimate 
the time-to-load, t^. When t = t^, x $ = x gm or x = 1, so that 



1 



A In 



cosh B(t.p + t ) 
t o 

cosh Bt Q 
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Under the further approximation (also applicable) that piston and 
slug damping forces are negligible, we have C->o and t Q +o, giving 



t f = 1, cosh -1 exp ^ (14) 

1/2 

|[1 - (k } + k 6 + k 9 )] (k 3 H 3 ) | 

where B 1 = 

2k 2 (l + H ] ) 



Each of the terms in this expression contain design parameters. In 
particular, it can be immediately seen that 



1/2 



2 m s D c \ 


\ 1 + H ! 


\ P /c , 


' ! 


(t) 


l H 3 


r/ivf .!-] 

[l D J P rJ 


1 


r 



(15) 



where H-| and H 3 are defined with Eq. (9) and P^ is an effective 
back pressure due to static drag and back pressure given by 



F F 

P = -2- + —5. + p 
*d A. A K b 
l c 

The value of B't^ given by Eq. (14) is plotted in Fig. 1. These 
expressions will receive additional attention in the discussion of 
the experiments. 
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t f [ a 

L ; 



^k 



2k 2 (1 + H X ) 



0 



1.5 



1 + h 2 

1 + H x 



Figure 2 . 
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2. Other closed form solutions : Other less realistic simplifica- 

tions can be made to Eq. (11) and these are listed here for completeness. 
Further details are given in Appendix B. 
a. Inertia effects only: Here, 

k 2 X (1 + H 1 + H 2 x) = a 
where a = 1 - (k-j + kg + kg) = 




The solution, developed in Appendix B, is: 

H, 



1/2 



t = 






2k 2 (1 + H ] ) / 1 + 



D(u) 



where D(u) is Dawson's integral [21] 



D (u ) = e 



“ 7 > : 



dt 



H, 



and, u = [fn (1 + x)] 



1/2 



Setting x = 1 we have, for the ram time 



[ 2k 2 (l + H ] ) 



t 1/2 



0| + H 



= 1 + 



1 



1/2 



1 + H l/\ H 2 



D(u) (16) 



This result is plotted in Fig. 2 for representative values of 
Ho 

1 + 



1 + H ] * 
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Note that is a measure of the nonlinearity contributed by the 
difference in acceleration of the liquid between injector and piston 
(H^ = 0 if D c = D^). Figure 2 shows the very small effect of this 
factor for the acceleration only case. Therefore, if this case is 
a reasonable approximation to Eq. (11), it is also reasonable to 
write 



l <2 x(l + H-j ) = a, so that 



t = 



• 2 k 2 (l 



1/2 

IV. 



and t. 



2k 2 (l + H 






1/2 



= 1 



(17) 



as seen in Fig. 2 . 

b. Fluid friction only: In this case, Eq. (11) becomes 



^3 (H 3 + H^x) = a 




so that 



(Appendix B) 




(18) 



Since is essentially equal to f c> the friction coefficient for 
the chamber, it is not likely to have a value larger than about 0.05. 

is of the order of magnitude of 1 (in the experimental apparatus 
it is about 3.3) so that it is unlikely that H^/H^ would be larger 
than 0.05 in which case t^ (a/k^H^)^ = 1.0124. Thus, it is seen 
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again that the nonlinearity contributed by the const. (x) term, as 
in the acceleration case, is relatively insignificant. With this 
approximation, we may write 

2 



^3 H 3 * = a, so that 



t = 




and t^ 




(19) 



This result is all the more useful in comparing the linearized 
versions of the acceleration only and fluid friction only cases. 

Thus , 

1/2 

r _ ram time, acceleration only 
ram time, fluid friction only 



1 i near i zed 



■ 



v I 



“1 / 






(20) 



which is a valuable ratio in analyzing the extent to which a system 
design is limited by inertia or by fluid friction effects. 

c. Completely linearized system: If the fluid friction is com- 

pletely neglected, along with H^, the result is the familiar linear 
damped spring-mass expression: 



k 2 x (1 + H-| ) + (k 4 + k 7 ) x + (kg + k g ) x 



with the solution 



-CO) t 

_ ^ n 



x = 



k 5 + k 8 



1 + 



sin r u _/T 

r. 2 n 

l - c 



t + a] 1 



( 21 ) 



where w 



2 k 5 + k 8 



. -1 /T 2 

a = tan -/ 1 - g 
_ k 4 + k 7 



n k 2 (1 + H ] ) and 2cto n k 2 (l+H, ) 
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Since no accurate values for kg, kg, k^, and k^ were available at 
the time of this writing (in the experimental apparatus, these values 
are essentially zero), this solution has not been analyzed in any 
great detail. In fact, our results indicate that these damping and 
spring terms do not significantly affect the performance of an LP6 
loading system when compared with the similar effects present in the 
compressible (gas) side of the system. 

In the following section of this report these closed-form solu- 
tions are compared with the experimental results. As has been men- 
tioned, the complete solution to the governing differential equations 
has been obtained by analog computer methods and is discussed in 
Section IV of this report. 
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III. EXPERIMENTAL PROGRAM 



A. DESIGN OF THE EXPERIMENTAL APPARATUS 

After a review of the available background literature and current 
research efforts, work was begun to design and fabricate a laboratory 
model capable of simulating the dynamic behavior of a liquid propel- 
lant feed system for a 20mm gun firing at a rate of up to 1,000 
rounds per minute. The selection of a 20mm bore diameter was dic- 
tated by considerations of laboratory capabilities and was not due 
to a desire to accurately simulate an operational system. The concept 
criteria coupled with the unique requirements for an accurate labora- 
tory simulation resulted in the establishment of the following initial 
design parameters: 



chamber volume (variable) 50-70cc 

maximum operating frequency 17hz 

water flow rate (maximum) 75.5gpm 

maximum chamber pressure 1500psi 

desired maximum time to fill chamber 15ms 



Due to the time constraints involved, it was felt that off-the-shelf 
items should be used whenever possible. The 20mm bore-diameter chamber 
was loaded with a brass slug weighing 93 grams which rode on two graph- 
ite filled Teflon sealing rings. The brass slug, which simulated the 
projectile, was cycled from the breech end of the chamber to the barrel 
end and returned to the breech end, completing one hypothetical firing 
cycle. To facilitate visual inspection of the bore and slug, the test 
chamber itself was fabricated from a three-inch O.D. Lucite cylinder, 

18 inches long, bored to a 20mm inside diameter and fitted with alumi- 
num end caps. The entire chamber assembly is held together with four 
tie rods and cap screws, as shown in Figs. 3 and 4. 
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Figure 3. Experimental Apparatus 




Figure 4. Gas Supply System 
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Because of the desire to vary the charge to mass ratio, a variable 
chamber volume was necessary. To accomplish this with one chamber, a 
volume control retaining rod was designed into the system. This brass 
rod, bored to allow gas to pass its length, was threaded through a 
plate which was attached to the barrel end cap holding the rod in the 
chamber. The rod, which has a Teflon disc on the end, not only estab- 
lished the volume of the test chamber, but provided a buffer stop for 
the slug at the end of its forward motion. Another Teflon buffer was 
affixed to the breech end cap to cushion the slug in return motion. 

The initial design concept started a simulated firing cycle with a 
slug at the breech end of the empty chamber, as shown in Fig. 5. This 
is the ready-to-ram position. The simulated propellant was then intro- 
duced, ramming the slug to the opposite (barrel ) end of the chamber as 
the chamber is filled. This was accomplished by applying gas pressure 
to a power piston which drove an injection piston. The injector piston 
forced the propellant past a flow check valve and into the chamber. 

This placed the slug in the ready-to-f ire or in-battery position. In 
an actual gun, the propellant would be ignited at this time in the cycle. 
Due to the laboratory constraints, the water expulsion system was used. 

With the slug at the in-battery position, ram gas pressure on the 
power piston was vented and return gas pressure applied to both the 
power piston and the muzzle end of the test chamber (Fig. 6). At the 
same time, gas pressure was applied to a pilot-operated check valve, 
the dump valve, to open it. With the dump valve open and reverse flow 
in the feed line eliminated by the flow check valve, the gas pressure 
on the muzzle end of the chamber forced the slug to the breech end. 
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Power Piston Injector P — i Test Chamber 
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Figure 5. Pneumatic control system in ram position. 
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Figure 6. Pneumatic control system in return position. 




expelling the propellant through the dump valve. At the same time, 
with the power piston and injector piston moving to the rear, another 
shot of fluid was drawn into the injector from the supply line. To 
repeat the ram stroke, the gas pressure was vented to atmosphere. 

Since the timing and synchronization of the gas and liquid flows 
were critical to the success of the experiment, a single control mech- 
anism was desired. With this in mind, an integrated gas and liquid 
system was designed. It was initially thought that balanced spool 
valves could be used to control one or both ends of the system. How- 
ever, obtaining off-the-shelf stock for such specialized requirements 
(the necessary low response times and high flow rates) presented sched- 
uling difficulties which were unacceptable. It was discovered that 
normally open and normally closed solenoid operated gas spool valves 
were "readily" available for operating pressures up to and including 
150psi and could be obtained with modifications allowing operation up 
to 250psi. 

A system was designed using one valve to control the ram gas pres- 
sure and one valve to control the return gas pressure. By using two 
valves, one normally open and one normally closed, gas pressure could 
be applied to either the ram or the return line with the other line be- 
ing vented by energizing or de-energizing the appropriate solenoids by 
means of a single relay. By connecting the pilot line of the dump valve 
to the return gas line, the entire feed system could be controlled with 
one electrical switch, or each valve could be energized and de-energized 
separately. This control system would also allow maximum flexibility 
for future efforts with little modification. Through an easily constructed 
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electrical timing circuit, automatic ram and return cycling could be 
accomplished for simulation of burst firing. 

B. CONDUCT OF THE EXPERIMENT 

The assembled experimental LPG simulator is shown in Fig. 7. A 
linearly variable differential transformer (LVDT) was manufactured and 
mounted next to the injector piston. This LVDT was attached to the 
connecting rod, between the power piston and the injector piston. 

The volume of the liquid being placed under pressure during each shot 
was measured by filling the system in the ready-to-ram position and 
then draining it into a graduated beaker. By measuring the displace- 
ment of the injector piston head with the LVDT, the volumetric rate of 
fluid injection during the ram stroke was obtained. 

Five pressure taps were drilled in the test chamber. Two of these 
taps were located as close to the breech end as possible, one at 20 
degrees from top center and the other at 20 degrees from bottom center 
in a counterclockwise direction, as seen from the breech end of the 
chamber. The remaining three taps were placed two inches, three and 
one-half inches, and five inches from the breech end, in line with the 
first bottom tap. 

It was discovered during operation that when the slug was returned 
to the ready-to-ram position, it would wipe some of the liquid out of 
these last three taps. In the next ram cycle, bubbles of gas which had 
replaced the liquid would distort the pressure response of these taps. 
To expedite the investigation, these three taps were plugged and only 
the taps closest to the breech were used. 
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Figure?. Chamber Disassembled 




Figure 8. Chamber Assembled 
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A 4-channel Hewlett-Packard 3960 Magnetic Tape Recorder was used 
to FM record desired data during system operation. Kaman diaphragm 
type (lOOOpsi) pressure transducers were connected to the two breech 
pressure taps on the test chamber. The pressure signals were processed 
with a Kaman Digi-Vit Readout Unit which also provided a visual 
(digital) display. Since the maximum output of the Digi-Vit unit 
was one volt, the breech pressure signals were amplified through 
Sanborn 8875A Differential Amplifiers prior to being recorded. The 
LVDT excitation signal was obtained by using a Hewlett-Packard Model 
200CD Oscillator. The output of the LVDT was passed through a voltage 
divider circuit and then recorded on the tape. 

A Brush Recorder (Mark 280) was used to obtain a visual display 
of the recorded data. By transcribing the desired signals on the 
Magnetic Tape Recorder at a tape speed of 15 feet per minute and 
playing them back into the Brush Recorder at 3-3/4 feet per minute, 
the time scale of the output was expanded by a factor of four on the 
Brush recordings (viz., from a real-time maximum of 200mm/sec to a 
delayed time maximum of 800mm/sec.). Two pressures and the LVDT 
reading were recorded for each shot. 

Before each shot at a new ram gas pressure, calibration points were 
marked on the Brush chart to establish the full-scale deflection pres- 
sures for each transducer during the shot. The LVDT full-scale deflec- 
tion was maintained at 40 divisions by adjusting the amplitude and 
frequency of the oscillator. A slug stroke of five inches was held 
constant throughout the experiment. 
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To eliminate the possibility of particulate contamination of the 
water, a 25 gallon tank of distilled water was connected to the dis- 
charge of the dump valve. This return line was elevated to place a 
hydrostatic head of approximately one psi on the discharge side of the 
dump valve to prevent gravity draining and to eliminate air ingestion 
in the chamber and connecting line when the dump valve was open. 

Data were obtained for nominal supply pressures ranging from 
50psig to the system limit of 220psig as indicated at the gas supply 
system regulator (Fig. 8). Further detailed descriptions of the ex- 
perimental apparatus and procedures are documented in a previous 
report [22]. 

C. EXPERIMENTAL RESULTS 

1. Chamber pressure : Ram displacement, chamber pressure, and ram 

gas pressure (driving pressure) were recorded for driving pressures 
between 50psig and 220psig in lOpsi increments. Typical results are 
presented in Figs. 9-12. In these figures, the recorder outputs have 
been arranged so that a vertical line crosses each trace at a corres- 
ponding time during the event. Event time increases from left to right 
with one centimeter equating to 12.5 milliseconds. Ram time has been 
defined as the time interval which starts when the ram gas pressure 
trace deviated from the horizontal to the point in time when any dis- 
cernable volumetric increase ceases, as shown in Fig. 9. 

After application of driving pressure to the power piston, a time 
lag was evident before the piston was set in motion. This time lag was 
largely due to static friction forces acting on the pistons and the slug. 
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Figure 9. Typical responses, 1^=50. 
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Figure 11. Typical responses, P r = 140. 
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Figure 12. Typical responses, P r = 210. 
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After an initial rapid increase of the driving pressure, these fric- 
tional forces were overcome, and injection of the fluid and, hence, 
the ramming of the slug proceeded. The relatively gradual increase 
in driving pressure after piston motion started has been attributed 
to losses in the upstream valving and connecting lines and the resultant 
inability of the supply system to provide an adequate flow of gas to the 
power piston chamber. 

During the ram stroke, the chamber pressure exhibited some charac- 
teristics which were not completely predictable by the analytical model 
in its present form. After an initial increase, the chamber pressure 
tended to decrease as the slug preceded the fluid column down the 
chamber. To facilitate the discussion of these pressure fluctuations, 
one must first observe the behavior of the displacement trace. In ad- 
dition to being a measure of the volume of fluid injected, it is an 
indicator of fluid velocity and acceleration behavior. As motion starts 
the near-zero slope of the trace implies a low velocity which increases 
as expected. The rapid initial rate at which the slope increases im- 
plies a rapid increase in the acceleration. As the ram stroke proceeds, 
this rate of increase of the acceleration diminishes and the curve be- 
comes more linear. 

As demonstrated by the analysis, the pressure difference across the 
system is mainly influenced by inertia forces and viscous forces. The 
inertial forces induce a pressure reduction proportional to the accelera 
tion. Likewise, the viscous forces endure a pressure reduction propor- 
tional to the square of the velocity. Following the initiation of 
motion, as the level of acceleration of the fluid and the slug mass 
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decreases due to friction, the difference in pressure between driving 
pressure and chamber pressure necessary to provide the acceleration 
level decreases with an attendant increase in the chamber pressure. 
Initially, fluid velocities are small so that the pressure difference 
is due mainly to this large but decreasing acceleration. As the accel- 
eration decreases, the chamber pressure increases until the velocities 
in the system become significant from a viscous friction point of view. 

As velocity increases, the loss in pressure due to these viscous effects 
overcomes the increasing pressure due to the decreasing level of accel- 
eration of the moving solid and fluid masses. The net effect at this 
point is a loss in chamber pressure. 

In several runs (Figs. 9 and 10, for example), the acceleration of 
the masses dies out after a few milliseconds. In these cases, the ex- 
pected continuous drop in chamber pressure, due to viscous effects, is 
observed after the initial pressure peak. In several other runs (Fig. 

11, for example), the acceleration effects are present throughout the 
run. In these cases (typical of higher driving pressures), the velocity 
of the moving masses increases through a large part of the event. Accel- 
eration levels decrease initially, as before, with the resulting initial 
rise in chamber pressure. Viscous losses are then evident with the 
corresponding decrease in chamber pressure. 

In a viscous fluid, pressure decreases are proportional to the square 
of the velocity while the sliding friction losses of the solid mass are 
more accurately modeled as being proportional to lesser power of the 
velocity. For a given change in pressure drop across the system, there- 
fore, there is a tendency for the slug to accelerate and decelerate more 
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rapidly than the fluid. In the present experiment, the slug apparently 
accelerates at decreasing levels throughout the ramming process; no 
deceleration of the slug was noted. Thus, the second and subsequent 
pressure decreases in the chamber, when they occur, are thought to be 
due to the inability of the fluid column to keep page with the slug. 

This motion appears to be oscillatory in the P c curves of Figs. 9-12. 

All of the effects discussed here are influencing the motion of the 
slug and fluid and these remarks are not meant to create an overly simple 
interpretation of the complex interactions involved. It is also impor- 
tant to note that the chamber pressure is measured at a fixed point at 
the breech end and that the amount of fluid upstream of this point de- 
creases with time. Thus, the length of the fluid column, upon which 
inertial and viscous forces are acting, decreases with respect to the 
location of pressure measurement. This change in length of the column 
is on the order of 30% in these experiments and tends to decrease pres- 
sure drops due to inertial and frictional effects as measured at the 
fixed pressure taps. The above discussion also assumes that the driving 
pressure is relatively constant for that portion of the event in which 
the forces mentioned are acting. The ram gas pressure did indeed behave 
approximately in this manner at lower pressures. As ram pressure was 
increased, however, the response deviated more and more from the assumed 
step input. 

2. Terminal behavior and ullage effects : The oscillatory be- 

havior of the chamber pressure after the slug came to rest was due to 
the expansion and contraction of the entrained gases in the fluid column 
[Figs. 9-12, for example). At the moderate driving pressures initially 
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used, the traces resembled second order damped responses and this 
resemblance prompted a further investigation of this portion of the 
ram cycle event. 

Viewing the dynamic portion of the system as picture in the sketch 

below, a force balance on the ram and injector pistons at the instant 

the slug reached the end of its stroke gives A (P -P ) = m x,- + 

p r c p i 

k fp*i + k sp x i ^ see Eq ‘ 




From the definition of the effective bulk modulus, , for a 
mixture of liquid and ullage gas, the flow due to compliance in the 
mixture may be written 



Vi ’ £ p c 



Neglecting the mechanical spring forces, k x.. 
expressions to eliminate x^. 



and combining these 



P = P 
c r 



~~ ("1 P „ + KpP ) 

2 „ pc fp c 



A P 6 e 



and, after Laplace transformation, and solution with P c (o) = 0 and 
P (o) = P Q , the pressure variation is found to be 
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= 1 + 



sin 



( 22 ) 





w n ^ 1 - x} t + cp 



where <J> = tan ^ - / 



2 





V 




and V- t = total fluid volume (liquid and ullage) 



Equation (22) for a damped linear system is analogous in form to the 
measured pressure variations following the ramming cycle. The natural 
frequency and damping coefficient can therefore be obtained from each 
driving pressure curve by measuring the appropriate amplitudes and the 
period of oscillation (see Appendix C). Using these experimental values 
of to n and r, , values for equivalent bulk modulus and friction factor were 
obtained. The variation in these calculated values were within the limits 
of uncertainty predicted for most cases [22]. The percentage of entrapped 
(ullage) gas in the fluid volume is then obtained by using this calculated 
value of Trapped gases were found to be on the order of 0.1 to 0.6% 
of the liquid volume during these tests. As the ram gas pressure was in- 
creased for successive runs, the oscillation of the chamber pressure be- 
came more erratic as can be seen by comparing Figs. 9 and 10 with Figs. 

11 and 12. The validity of the results obtained for u) n , and k^ 
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at these higher pressures is therefore suspect, even though the 
general form of the curve repeats at all driving pressures. 

In some of the events, the chamber pressure curve indicates a 
pressure below atmospheric immediately after the slug reaches the 
end of its travel. These initial low-pressure oscillations are 
tentatively attributed to slug/fluid interactions as the slug comes 
to rest against the buffer. Note that following the stopping of the 
slug, the initial pressure change is always a steep increase, indi- 
cating that the fluid column is initially compressed against the 
stationary slug. 

It was hoped that a value of equivalent bulk modulus could be 
obtained from the displacement trace to correlate the value obtained 
using Eq . (22 ) . However, electrical noise in the LVDT circuitry 
produced a signal with such little definition that the extremely 
small volume changes associated with trapped gas effects were not 
distinguishable. 

3. Ram time measurements: The data displays of ram-time versus 

ram gas pressure (Fig. 13) portray some rather interesting charac- 
teristics. The general trend and validity of the curve is confirmed 
by the fact that the data presented were taken on several occasions 
and repeatedly at selected conditions. It should also be noted that 
for driving pressures above the moderate range (p f 55 1 50ps i g ) the 
decreases in ram-time are small for increases in ram gas pressure. 

The scatter of the data can be attributed to the inherent uncertainty 
encountered in any experimental data taking system, coupled with pos- 
sible erratic behavior of the pneumatic equipment. Any uncontrolled 
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Figure 13. Ram Time vs. Ram Gas Pressure - Experimental 



variations in equivalent bulk modulus for the different tests did not 
appear to significantly influence the ram-time results. An effort was 
made throughout the experiment to keep the amount of entrapped gas small 
and relatively constant. The data points conforming to a consistent 
trend would indicate that this was accomplished. 

4. Comparison with predictions : The experiments discussed herein 

have been conducted with the following values of the design parameters: 



c 

C/ITI 



D. 



1 



l 



Wc 

k fp =k fs 

L. 

1 



m 

P 

m s 

n 

p d 

p 

r 

t 



sm 



2 

chamber cross-sectional area, in 


(0.49) 


charge-to-mass ratio 


(0.43) 


injector diameter, in. (=D r > in this case) 


(1.50) 


connecting line diameter, in. 


(0.57) 


chamber diameter, in. 


(0.79) 


flow friction factors 


(0.03) 


mechanical damping coefficient, lbsec/in. 


(.001) 


injector length, in. 


(1.60) 


connecting line length, in. 


(16.0) 


2 

mass of injector piston, lb sec /in. 


(5.2 x 10" 3 ) 


2 

mass of slug, lb sec /in. 


(5.3 x 10' 4 ) 


number of connecting lines 


(1) 


equivalent pressure drop due to drag, lb/ir 


o 

i. (empirical) 


2 

ram gas pressure, lb/ in. 


(control! ed) 


time, sec. 


(variable) 


ram-time, sec. 


(variable) 


slug position, in. 


(variable) 


maximum (final) slug position, in. 


(5.00) 
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The sudden decrease in diameter from injector to connecting line gives 
a minor loss coefficient (see Eq. (4)) of 




The entrance to the chamber is constricted by an orifice-type reduc- 
tion with contraction ratio (C c ) of 0.92. The coefficient for this 
is (with C d =0.61 ) 




and = 

These design factors lead to the following values of the relevant 
parameters for Eqs. (9) and (10): 



k 2 - 5.41 x 10“ 3 /P r sec 2 

k 3 = 7.36 x 10“ 3 /P r sec 2 

k 4 = 7.85 x 10“ 4 /P r 

k 5 = k g = 0 

k ? - 1.02 x 10“ 2 /P r 

H-j = 3.44 

H 2 = 0.39 

H 3 - 3.21 

H 4 = 0.03 

The value of the effective back pressure, P^ (see Eq. (15)) is obtained 

empirically as is discussed in subsequent paragraphs. 
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a. Pressure and displacement variations with time. 

Figures 9-12 show the results of the analog computer solu- 
tions to the analytical model for the conditions under consideration. 
The value of P^, the effective back pressure, was set at 46psi except 
for the P r = 50psi run (Fig. 9) where it was necessary to set 
P^ = 41 ps i . The logic of these selections is discussed in a subse- 
quent section. In obtaining the fit of the analytical model to the 
data, as was done in Figs. 9-12, it was necessary to partition P^ 
into portions acting across the injector piston and across the slug. 
That is, 



P d + + P b P di +P ds 



where P ..= and P 



di A. 



ds A, 



+ P L 



In the four cases illustrated, it was found that the static drag was 
almost entirely due to the injector piston. The semi-empirical re- 
sults are as follows: 
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This information is of obvious use in pointing out the area to seek 
reductions in P^. 

The agreement of the analysis with the experimental measurements 
shown in Figs. 9-12 is remarkably good when it is realized that (1) 
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the step input assumed in the model is only crudely approximated by 
the experimental system and, (2), the model does not include compres- 
sibility effects. The latter discrepancy is believed to be the reason 
that the analytical model does not accurately follow the observed 
chamber pressure oscillations. These oscillations, caused by the 
accordion-like motion of the liquid between the injector piston and 
the slug, require the accounting for compressibility effects in the 
liquid. A careful examination of these results has shown that the 
deviations from the analytical model are also present in the dis- 
placement trace, although this is not obvious in the figures. The 
effect of the non-constancy of driving pressure, P^,, is partially 
removed by the semi -empiricism present in the selection of P.^ and 
P s d> as discussed above. Finally, it should be noted that the cases 
presented in Figs. 9-12 are randomly selected from a large number of 
runs. The entire body of original data is available upon request, 
b. Ram-time variations with driving pressure. 

Equation (14) gives a predicted ram-time dependency based 
upon a number of simplifying assumptions. Repeating the expression 
here, 




and 



. - k 2 (1 tH l ) 
k 3 H 3 
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With the experimental conditions previously described, we have 



.156 
< P r - P d> 



1/2 



and 



with the result that 



A = 1.017 



i_ = 256 

■ f (p r - p d ) ,/2 



(msec) 



The agreement of this expression with the experimental data is shown 
in Fig. 14. To obtain this agreement, it has been necessary to account 
for the approximately 8 millisecond delay that is observed in the ex- 
periments. In addition, it is seen that a good estimate of P^, for 
the system used, is 46psi. This remarkable agreement indicates the 
essential correctness of the model when used within the limits imposed 
by the assumptions. The main uncertainty in the model is about ±10% 
in the value of H^. This uncertainty is equivalent to about ±2% in t^. 



5. General observations . The analytical model is seen to be well 
supported by the experimental data. It is therefore reasonable to 
utilize the analytical model to extrapolate the experimental results 
to cases of special interest. In particular, the analysis can be used 
to indicate the areas of redesign likely to be most profitable and to 
obtain preliminary predictions of LPG loading system performance under 
a variety of design conditions obtained through the specification of 
values of parameters different from those analyzed to date. These 
parameters include caliber size, charge- to-mass ratio, and projectile 
mass (or propellant density). Results of such parametric studies are 
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Figure 14. Ram time vs ram gas pressure; comparison with theory. 



reported in the following section. 

Finally, it should be pointed out that, as expected, the analytical 
model only partially predicts the instantaneous behavior of the liquid 
and solid masses during loading. If knowledge of the mean behavior 
with respect to time is insufficient information, it is imperative 
that the model be expanded to take the driving system (compressed gas) 
into account and to allow for compressibility effects in the liquid. 
Some further comments on this subject are offered at the conclusion 
of the report. 
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IV. DESIGN CONSIDERATIONS 



In this section the analytical model is utilized to demonstrate 
those design decisions that are most influential in determining the 
performance of LPG propellant loading systems. For the purposes of 
the following remarks, "performance" is taken to mean ram-time. Other 
criteria (such as chamber pressure levels and system volume) are also 
amenable to prediction by the analytical model but ram-time is con- 
sidered here due to its wide acceptance as a quantity to be minimized 
in any LPG system. 

A. REDESIGN OF THE LPG SIMULATOR 

Appendix D describes the analog computer program that has been 
developed to solve the governing equations (9) and (10). This program 
has been used to predict the behavior of the existing laboratory system 
and to point the way toward the most profitable areas for redesign em- 
phasis. 

Figures 15a-e show the computer outputs for slug displacement, 
velocity, acceleration, chamber pressure, and slug pressure, respec- 
tively, for a range of reservoir pressures from 100 to 250psig. The 
values of the parameters are those presented previously in Section 
III.C.4 of this report. The parameters k^ and k^ are mainly dependent 
upon operational specifications (slug mass, caliber, charge- to-mass 
ratio, etc.) and are therefore not suitable for redesign. The other 
k. parameters are negligibly small in the existing system. The 
parameters and H^, representing inertia and viscous effects, 
respectively, are representative of two central areas for redesign. 
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With H-j = 3.44 and Hg = 3.21, we see from Eq. (20) that R = 1.42 and 
the system is therefore inertia dominated from a ram-time point of 
view. The third parameter that might be improved is the static drag 
on the system, represented by the effective back pressure P^. 

Figures 16a-c show the effects of the parameters , Hg, and 
at a reservoir pressure of P^ = 150psig. Each of these parameters 
have been swept through four-fold reductions from the existing design. 
From these figures it is seen that ram-time decreases of 3.5, 5.5, 
and 2.5 milliseconds are predicted for four-fold decreases in P^, 

H-| , and Hg, respectively . It is apparent, therefore, that redesign 
priorities should be inertia reduction (H^), static drag reduction 
(P d ), and viscous loss reduction (H^ ) , in this order. Inspection of 
the parameters H-] and Hg will show that a reduction in connecting 
line length (L^) and/or an increased connecting line diameter (D^) 
will decrease both H-j and Hg. A significant reduction in is 
attainable through the decrease of piston mass (nip) or the increase 
of injector area (this term is already rather small, however, rela- 
tive to the third term in H-j). Hg can be significantly reduced 
mainly through shortening of line length, increasing the line diame- 
ter, and, most importantly, reduction of minor losses. In the 
current system, minor losses make up about 85% of Hg (thereby point- 
ing up a discrepancy in terminol ogy) . The effective back pressure 
is apparently due mainly to drag at the injector end of the system. 
This term is also influenced (due to inadequacies of modeling of the 
input pressure) by a non-constant reservoir pressure. 
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Figures 17a-e show predicted performance curves calculated by 
assuming a two-fold reduction in all of the redesign parameters. 

That is, in these calculations, = 23psig, H-j = 1.72, and = 1.60. 
Referring to the P^ = 150psig curve, we see a ram-time reduction of 
about 7.5 milliseconds. This appears to be about one millesecond 
more than that calculated by simply adding the separate reductions 
observed in Fig. 16 (thus, the non-1 inearity of the problem is mani- 
fested). Figure 17 shows that if the two-fold parameter reductions 
are attainable through redesign, a ram-time of 15 milliseconds should 
be attainable with a reservoir pressure of about 200psig. Allotting 
25% of the total firing cycle to liquid loading, this would permit a 
rate of fire 1000 rounds per minute. 

In concluding this section, it should be noted that all delays in 
loading, other than those due to actual motion of the injector, liquid, 
and slug, have been omitted from this discussion. The 8 millisecond 
static delay observed in our experiments would, of course, be unaccept- 
able in an operational system. With this caveat, however, redesign 
recommendations, in order of importance, are: 

1. Reduce connecting line length. 

2. Increase connecting line diameter. 

3. Decrease piston mass and static drag. 

4. Seek out and reduce minor losses. 

B. PRELIMINARY DESIGN OF A 5- INCH SYSTEM. 

In this section, the analytical model is exercised in order to 
demonstrate its utility in making first-cut design decisions applicable 
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1 Figure 17. Predicted performance, LPG redesign. 



62 




Figure 17 (concluded) . 
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to liquid loading systems for guns of different scale. Before pro- 
ceeding, it is noted, once and for all, that the validity of these 
predictions is vitally dependent on the extent to which the basic 
analytical model approximates the system under design. In particular, 
any current injector design will probably require a more complex 
mathematical description than that given by Eqs. (1) and (2). 

For the purposes of this discussion, we shall rely mainly upon 
Eq. (14) which describes the semi-linear case in which H^/Hg, H^/ ( 1 ) , 
and (k^+kg) are all negligibly small. These are desirable design con- 
straints and are relatively easily obtained. In any case, in sections 
11. A. 2. a. and b., the effects of departures from these conditions have 
been shown to be small. Repeating Eq. (14) in a slightly reworked form 



k ? (l + H, ) ? 

where A = ■ . „ — — = R /2 

k 3 H 3 

Note the physical interpretation of the parameter A: it is a measure 

of inertia effects relative to viscous effects. The function of A 
in Eq. (23) is shown in Fig. 18. 

It is now necessary to postulate a few design parameters descrip- 
tive of a 5-inch base! ine LPG propellant loading system. It is to 
be emphasized that this baseline design need only be approximate 
since it is the function of the model to point up profitable depart- 
ures from the baseline. For this example, the baseline system is 



k 3 H 3 



A?) 



i/^ 



A cosh'l exp (1/A) 



(23) 



r _ 
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Figure 18. 



described as follows: 



Chamber diameter 
Charge-to-mass ratio 
Propellant density 
Slug mass 
Driving pressure 



D c = 5.0 in. 
c/m s =0.5 

p = 1.3 x 10 ^lb sec^/in.^ 
m = 80 lb = 0.207 lb sec^/in. 
P r = 250 psig 



= 
m 

s 



= 1 . 0 , 



= °- 5 > 



= 0 . 1 , 



r 5 -°- 



n = 2 

57 " '•». 
1 



xt* K0 - 

sm 



'-p 

4 = '- 0, 

sm 



£ k n,i - 20 



The length of slug travel can be calculated from 

m 



sm 



s(c/lTl ) /in c ■ 
s' = 40.6 in. 



pA. 
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The pertinent parameters are found to be 

k 2 = 1.71 x 10" 3 sec 2 
k 3 = 3.47 x 10~ 3 sec 2 
H ] = .0625 + .125 + 6.25 = 6.44 
H 3 = .0009 + 23.44 + 2.46 = 25.90 

Note the dominance of the connecting line contributions, 6.25 and 23.44, 

to H-| and H^, respectively. From these values, we calculate 

A = 0.141 and R = 0.53 

so that this design is seen to be dominated by viscous effects. From 
Fig. 18, f(A) = 1.1 and the ram time is t^ = 0.35 sec. With the criter- 
ion of 25% of the total cycle time available for liquid loading, this 
baseline ram-time will allow a rate-of-fire of about 43 rounds-per-minute. 

There is a limitless number of perturbations on the baseline design 
that can be performed in order to examine their effect upon rate-of-fire. 
Judicious choice of these perturbations awaits the specification of more 
exact performance and design criteria. As examples of what can be done, 
however, it is noted that this particular baseline design is severely 
constrained by inertia and viscous losses in the connecting lines be- 
tween injector and breech. Thus, if the number of 1-inch connecting 
lines is doubled (to n = 4), we calculate H-| = 3.31, and = 8.32, so 
that A = 0.255 and R = 0.714. The resulting maximum rate-of-fire is 
about 71 rounds-per-minute, a 65% improvement. As another example, if 
the injector diameter is constrained so that D /D. = 1.0, the resulting 

C 1 

rate-of-fire, for the n = 4 case, is about 55 rounds-per-minute, a 29% 
decrease. 

As a final demonstration of the analytical capability thus far dev- 
eloped, the baseline 5-inch design has been modeled on the analog compu- 
ter for various driving pressure levels. The results are shown in Figs. 
19a-e. 
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Figure 19. Predicted performance of 5-inch baseline design. 
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V. SUMMARY 



A. CONCLUSIONS 

An analytical model has been developed to predict the time-dependent 
behavior of the moving liquid and solid masses during an LPG loading 
cycle. This model has been tested against experimental results obtained 
with an LPG laboratory simulator and has been found to be adequate for 
the prediction of system performance. The model is specifically con- 
strained in its use by the assumption of an incompressible liquid and 
a uniform step input in driving pressure. 

An experimental LPG simulator has been designed and fabricated to 
assist in verifying analytical predictions and to identify design and 
performance phenomena requiring particular attention in future LPG 
developments. Measurements made with the LPG simulator, together with 
guidance provided by the analytical model, have provided valuable in- 
sights into the behavior to be expected from a given LPG feed system 
design. Specific unanswered questions relate to the effects of large 
trapped gas volumes and the existence and probable effects of cavita- 
tion in the liquid slug during the loading cycle. 

The analytical model has been shown to be useful in the prediction 
of the influence of design specifications such as caliber, projectile 
mass , charge- to-mass ratio, propellant properties, and system geometry 
and scale. It is anticipated that the analysis, further substantiated 
by experimental observation and improved as required for future appli- 
cation, will provide a valuable tool in the design of prototype and 
operational LPG systems. 
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B. RECOMMENDATIONS 



1 . Experimental 

a. Incorporate design improvements as suggested by the analyti- 
cal model. Test the modified LPG simulator to verify the 
expected performance increases. 

b. Conduct tests to verify other parametric dependencies pre- 
dicted by the model. Those most easily evaluated include 
propellant properties (simulate a liquid propellant), 
charge-to-mass ratio, and slug mass. 

c. Investigate the effects of large volumes (up to 10%) of 
ullage gas. 

d. Determine the existence of cavitation and, if it is found, 
evaluate the effect of cavitation upon system performance. 

e. Develop the simulator to the point that it can be operated 
in the fully automatic mode. Conduct tests appropriate 
for the evaluation of the system in this mode. 

2. Analytical 

a. Continue to refine the model. The difficulty of including 
compressibility effects may preclude this improvement un- 
less future requirements justify the necessary effort. 

b. Program the governing equations for solution on the digital 
computer. This will be necessary if the model is extended. 

c. Model upstream driving systems in order to obtain accurate 
inputs to the existing model. This will be necessary if 
the analysis is to provide performance predictions for the 
fully automatic mode of operation. 

d. Monitor LPG program developments with a view to subjecting 
emerging designs to evaluation by the model. 
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APPENDIX A 



a 

A 

C 



(c/m $ ) 

D 

f 

F 






m 

k] • • kg 



m 

n 



P 

P 

P 

R 

t 

t 



b 

d 

f 



NOMENCLATURE 

parameter defined in conjunction with Eq. (16). 

cross-sectional area 

contraction coefficient 

discharge coefficient 

charge-to-mass ratio 

diameter 

fluid friction factor 
static drag force 

parameters defined in conjunction with Eqs. (9) and (10) 
damping force coefficient 
spring force coefficient 
minor loss coefficient 

parameters defined in conjunction with Eqs. (9) and (10) 

length 

mass 

number of connecting lines 
pressure 

back pressure on slug 

effective back pressure defined in conjunction with Eq. (15) 

ratio defined in Eq. (20) 

time 

final time, ram-time 
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x displacement, position 

x sm final slug displacement 

8 effective bulk modulus 

e 

? damping factor 

p liquid density 

w n natural frequency, radians/second 



Subscripts 

c refers to chamber 

i refers to low pressure (liquid) side of injector; also 

counting index 

l refers to connecting line(s) 

p refers to injector piston 

r refers to high pressure (reservoir) side of injector 

s refers to slug 



74 



APPENDIX B 



CLOSED FORM SOLUTIONS TO THE GOVERNING 
DIFFERENTIAL EQUATION 



We seek solutions to the equation 

■2 



k^x ( 1 + H-j + Hg x) + k^x (H^ + H^x) + (k^+k^)x + (Xg + kg)x = a (Bl) 



x(o) = x(o) = 0 



1. Semi -linear 



h 4 h. 

Let — << 1 , jq qj — << 1 , kg + kg << 1 , so that 



kg (1 + H-j ) x + kgHgX^ + (k^ + k^) x = a 



Let z = x, then kg(l + H-j ) z + k^H^z 2 + (k^ + k 7 ) z = a, and 
dz dt 



2 " k 2 (l + H-j ) 



a - (k^ + Y-j)z - kgHgZ 



and, after integration and application of the limits 



q 



tanh" 



•2k 3 H 3 z - (k 4 + k 7 )' 






tanh 



-1 -(k 4 +k 7 ) 



f 



“ k 2 (l + H-j ) 



where q = (k^ + k 7 )^ + 4akgH 3 
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Upon rearrangement and solution for z = x, we have 



x = AB tanh B(t + t ) - C (13) 



where the constants are defined in the main body of the report. 
Expressions for x(t) and x(t) are obtained through integration 
and differentiation of this expression. 

2. Inertia effects only . Equation (B1 ) becomes 

l <2 x (1 + H.| + H^) = a 

Multiplying by 2x 



2x x = 



2xa 



k 2 (l + H 1 + H 2 x) 



(B2) 



or 



so that 



./•2x _ 2a dx 

d(x ' k 2 (1 + + Hgx) 

x - + rr« 

Ho " 1/2 



1 



1 V2 



Let u = [£n (1 + y - ^ x)] so that dx = ^ 



2(1 + H } ) u 2 



u e du 



(B3) 



Equation (B3) becomes 



2(1 + H-j ) u 2 



e du = 






2a 



k £ H 2 



dt 



or 



v? [^’] //.. 



= t 



and, since e 



-u 



H, 



upon rearranging. 



1 + 



1 + H 



1 



s ( 1 + tVh, x ) . 



2k 2 (l + 




1 + H 
~HZ~ 



1/2 






D(u) 
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-u 2 f U 2 

where D(u) = e u J e z dz is Dawson's integral [21]. 
o 

3. Fluid friction only . In this case, Eq. (B1 ) becomes 

k^x 2 (H^ + H 4 x) = a 

(H 3 + H 4 x ) 1/2 dx = (^) 1/2 dt 

and, upon integration and evaluation of the limits, 




Note that since we have reduced the order of the basic equation, the 
initial condition x(o) = o cannot be satisfied in this formulation. 
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APPENDIX C 



NATURAL FREQUENCY AND DAMPING 

The natural frequency, w , and the damping coefficient, c , 
can be calculated by measuring the quantities indicated on the sketch 
be! ow: 




C = 



1 ln[a/A] 



+ 1 



CO = 

n 



2tt 



V 



T V 1 - ? 



where to is in rad/sec and T is in seconds, 
n 



A. Effective Bulk Modulus. 



to„ 



B e = 



a Al 

IlL V + 
p t 

m V. 9 
p t 2 



<o„ 



A. 
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GO 

n 


natural frequency, rad/sec 


(variable) 


a p 


. 2 

piston area, in 


(1.767) 


% 


2 

mass of piston, .^ ec 


(5.2 x 10" 


v t 


3 

total volume of fluid, in 


(21.389) 




equivalent bulk modulus, psi 


(variable) 



The values in parentheses are for this experimental setup. 



) 



3„ = 0.0354 .2 
e n 



B. 



Friction Factor, 



k 



fp 



c = J ■ 

2A p \ 

k fp = 2% “n c 

k- = 0.01035 co c 

f p n 



C. Trapped Gas (Ullage). 



1 _ 

e e 



h 6 g V t 



3 g equivalent bulk modulus 

3^ liquid bulk modulus 

3^ gas bulk modulus 

V„ total volume of gas in fluid 

9 



For a perfect gas 



where y ratio of specific heats 

P absolute pressure 
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For air at moderate to low pressures, it can be assumed that air 



approximates a perfect gas 




For water 

^ « 3 x lO^psi 



Is. = 1.4P / 3 x 10^ - 6 e V 1004 
V t \3xl0 5 6 e / 

As previously reported in [22], the trapped gas volume thus calculated 
ranges from 0.1% to 0.5% of the total contained volume V^. 
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APPENDIX D. 



SOLUTION OF THE GENERAL 
GOVERNING EQUATIONS. 



The governing equations, developed in Section II of this report, 
are repeated here: 

1 - P $ = k-| + k^ x(H-| + H^x) + kgX 2 (Hg + H^x) + k^ x + kgX + kg (D-l ) 

P $ = k 2 x + k 7 x + kg x + k g (D-2) 

The solution of these equations, on the analog computer, is developed 
in this Appendix. 

The variables in Eqs. (D1 ) and (D2) are non-dimensional in all 

• 2 

respects except for time. For instance, the dimension of x as it 

-2 

appears in Eq. (D1 ) is sec . Non-dimensional ization with respect 

to time is now advantageous for scaling purposes. 

Letting non-dimensional time r = gt and defining 3 such that 

B 2 k 3 = 1 (D-3) 

we have, after substitution, elimination of P , and rearrangement 



*■&[ 



1 - 



(k-j + k 6 + x ^ H 1 + H 2 ” *^3 + H 4 x ) " 

3 



6 (k 4 + k 7 ) * ] ' if < k 5 + k 8 )x 



(D-4) 



Note that although the notation has not been changed, all variables 
in Eq. (D-4) are completely non-dimensional. The expression for 
the acceleration has been written in the form of Eq. (D-4) in order 
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to illustrate the scheme of programming for the analog computer. 



An auxiliary relation for the slug pressure P $ is 

k ? 

P s = kT * + k 7 6 * + k 8 x + k 9 



(D-5 > 



The chamber pressure, P c , can be developed by adding to P s the 
pressure drops due to inertia and viscous effects in the chamber. 
The non-dimensional result is 



Equations (D-4) through (D-6) are shown in schematic form in 
Figs. D-l(a) and D-l(b). Figure Dl(c) illustrates a modification 
necessary when the design is dominated by viscous effects (^>10) 
such as in the 5-inch baseline example of Section IV. It should 
be remarked that the logic of Fig. D-l is based upon unit program- 
ming (reference voltage = 1 unit) and that the schematic illustrates 
the manner in which the equations were programmed for solution on 
two full -expanded EAI TR-20 10-volt analog computers. The primed 
component numbers in Fig. D-l(b) refer to the slaved computer and 
numbers without primes refer to components of the master unit. The 
table following Fig. D-l lists the parameters determining the 
potentiometer settings in the analog program. 




where 




(D-6) 
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Figure Dl. Wiring schematic for analog solution. 
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Figure Dl (concluded) . 



Pot Number 


Parameter 


1 


(k 4 + k 7 )e 


2 


(k g + kg) Ckg/kg) 


3 


H 4 /10 


4 


Hg/l 0 


5 


h 2 /10 


6 


H-j/10 


7 


k 2 /k 3 


8 


kg/lOkg 


9 


k ? 6/10 


10 


k 8 


11 


k 2 /k 3 


12 


k 9 


13 


1 _ (k-j + kg + kg) 


15' 


(k 2 /k 3 )(c/m s ) 


17' 


f 




s 



Table D-l. POT SETTINGS 
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